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Abstract 

Some particular examples of classical and quantum systems on the lattice are 
solved with the help of orthogonal polynomials and its connection to continuous 
models are explored. 
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1 Introduction 

Recently there is an increasing interest on discrete models in classical and quantum 
physics. Numerical calculations use computational technics to solve differential equa- 
tions [2]. Lattice field theories have become a very powerfull tool to avoid infinities in 
perturbative methods, and to obtain exact solutions of the field equations [7]. Ponzano- 
Regge calculus introduced in gravitational field is equivalent to Penrose spin network that 
discretizes riemannian manifolds [1]. Statistical mechanics has been working from the 
begining with lattice approximation. 

The orthogonal polynomials of discrete variable [8] offers a new approach to these 
models. They are exact solutions of difference equations from which raising and lowering 
operators, eigenvalue equation, symmetries and constant of motion can be calculated. If 
a physical problem is not given in a discrete form, this can be guessed if we know the 
differential equation of some orthogonal polynomial which is the continuous limit of the 
diference equation. 

We present some simple examples of discrete models in classical and quantum systems, 
that can be solved by the above method. In the last case, special care has to be taken 
in the construction of Hilbert space when the basis are the orthogonal polynomials of 
discrete variables. 

2 Orthogonal polynomials of discrete variable 

A polynomial of hypergeometric type Pn{x) of discrete variable x satisfies two fundamental 
relations from which one derives raising and lowering operators [6]: 

i) Difference equation: 

a{x)AVPn{x) + r{x)APn{x) + A„P„(x) =0 (1) 
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where a{x) and r(x) are polynomials of, at most, second and first degree, respectively, 
A(V) are the forward (backward) operators and A„ is the eigenvalue corresponding to the 
eigenfunction Pn{x). 

ii) Three term recurrence relations: 

(^) (2) 

iii) Raising operator 

a{x)VPn{x) 

where 

r„{x) = 

iv) Lowering operator 



= ^rn{x)Pn{x) - ^Pn+l{x) 

nr' In 



(3) 



t{x + n) + aix + n) — o{x) 



[U{X) ^ t{x)) ^Pn{x) 

2n ^ 



n A2n+i 2n 



(4) 



These polynomials are orthogonal with respect to the weight function p{x\ particular 
examples of which are the Kravchuk, Meixner, Charlier, and Hahn polynomials. 

Prom the orthogonal polynomials of discrete variable we can construct the correspond- 
ing orthonormal functions of discrete variable, by the definition 



= d^^ \J p{x)Pn[x) 

where dn is some normalization constant. They satisfy 
i) Difference equation 



(5) 



\/{u{x) + r(,T)) a{x + \)^n{x + 1) + -J{(y{x - 1) + t{x - 1)) a{x)^rix - 1) 

-(2a(x) + r(x))0„(x) + A„</.„(x)-O (6) 



ii) Three term recurrence relation: 

— an—^(t>n+l[x) + ^7n-^0n-l W + l^\Pn " x)(l)n[x) = 

iii) Raising operator 
"A„T„(x) 



(7) 



n r' 



^n{x) + Vi'^i^ - 1) + ^(^ - 1)) '^{x)Mx - 1) 



A2n , / X 

= -^Oin—J—<Pn+l{X) 

In dn 



(8) 
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iv) Lowering operator 

L~{x, n) = 



KrJx) , A 



'r^ + K+ '-^{x - Pn) - cr{x) - t{x) 



n tL 2n 



+ ^{a{x)+T{x))a{x+l)(l)n{x + 1) = (9) 

The difference equation i) written in the form H{x,n)(f)n{x) — 0, can be factorized 
with the help of the raising and lowering operators: 

L~{x, n + l)L^{x, n) = fi(n) + u{x + 1, n)H{x, n) (10) 
L"'"(a;, n)L~{x, n + 1) = ii{n) + u{x, n — l)H{x, n + 1) (11) 

where /x(n) = — - — — - an7n+i, u[x, n) = a{x) 

2n 2n + 2 nr' 
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3 The quantum harmonic oscillator of discrete 
variable 

We start from the orthogonal polynomials of a discrete variable, the Kravchuk polynomials 
kn\x) and the corresponding orthonormal Kravchuk functions [5] 



Kif>\x) = d;,'.M^k^^\x), (12) 

where rf^ = :^T(^~^(w)" is a normalization constant, p{x) = ^^^^j^Zx^ivo)'^ is the weight 
function, with p > 0, > 0, p + q— 1, x = 0, 1, . . . — 1. 
From the difference equation 



^Jpq{N ~x){x+l)K^;^{x + 1) 

+ ^Jpq{N - x + l)xK^^^ {x-l) + [x{p -q)-Np^n\ KjP^ (x) = 0, ( 13) 

and the recurrence relation 



y/pq{N -n){x + l)K.^ 



(p) 

n+l 



X 



+ ^/pq{N -n + l)nK'^J^l^{x) + [n{q - p) + Np - x] K^\x) = 0, (14) 
we construct raising and lowering operators 

L+{x, n)Klf\x) = pq{x + n- N)Klf\x) 



+ ^Jpq{N -x+ l)xK^^^ - 1) = \/pq{N -n){n + l)K^l^ (x) , (15) 

L-{x, n)KjP\x) = pq{x + n- N)K^\x) 

+ y/pq{N -x){x + l)K^\x + 1) = ^Jpq{N -n + l)nK^^\{x). (16) 

As in the general case, these operators factorize the difference equation 

L+(a;,n - l)L~{x,n) ^ pq{N - n + l)n + pq{x + n - 1 - N)H{x,n), (17) 

L~(x, n + l)L'^{x, n) = pq(N -n)(n + 1) +pq(x + n + l- N)H(x, n). (18) 



In order to justify the name of quantum oscillator of discrete variable we substitute 
X — Np + ^j2Npqs, and take the limit ^ oo in the former expressions. We get 

where Hn{s) is the Hermite polynomial of continuous variable. The difference and re- 
currence relations for the Kravchuk functions becomes in the limit the differential and 
recurrence relations for the normalized Hermite functions ijj{s) which are the solution of 
the quantum harmonic oscillator. 

In the limit the raising and lowering operators becomes the creation and annihilation 
operators 

L-M KtH.) ^-.^ i {. - - (20) 

-±= L-(.,n) ir<^)(x) ^-.^ i {. + ^ am (21) 

The commutator of the raising and lowering operators are closed under the SO{i) algebra. 

^ [L, L+] K^\x) = (l - y) K^\^) ^ L'K^\^) (22) 

which in the limit becomes [a, a+] ip{s) = ip{s) 

Similarly the anticommutation relation of the raising and lowering operators becomes 
in the limit the hamiltonian of the quantum oscillator 

^ [L, L+] Klr\x) = J + 1) - (j - nr) K^\x) 

— y (aa+ + a+a) ^Pn{s) = (2n + l)^/'n(s) (23) 

4 The hidrogen atom of discrete variable 

We define the orthonormal Meixner functions [5] 



M(^'^) {x) ^ d-' ^r^) mi^''^) {x) (24) 

where rrin''^^ (x) are the Meixner polynomials, 

n!r(n + 7) , , /i^r(,T + 7 + 1) 



t,-{i-t,)^T{^y ' r(a; + l)r(7) ' 



and 7, are real constants 0</i<l, 7>0. 
They satisfy the orthogonality condition 



j;il4-)(x)il^7-)(x)-^ = 5„., 



and the following properties: 
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i) Difference equation 



' IJ,{X + 7)(,T + + 7) 



X + 7 + 1 



+ VM^T7)^M„(x - 1) - [ii{x + 7) + X - n(l - fj,)] Mn{x) = 
ii) Recurrence relation 



(25) 



- V/^(n + 7)(n + l)M„+i(a;) - ^/"(n + 7 - l)nM„_i(x)+ 
+ (/ix + /xn + /X7 + n — x) M„(x) = 

iii) Raising operator 

L+(x, n) M„(x) = + 7 + n) M„(x) + + 7)x M„(x - 1) 

= V//(n + 7)(n + l) M„+i(x) 

iv) Lowering operator 

L~{x, n) Mn{x) = -ii{x + 7 + n) M„(x)+ 



(26) 



(27) 



+ 



y X + 7 + i 



(28) 



Notice that we have omitted, for the sake of brevity, the superindices (7,/u) in M„(x). 

In order to make connection between the Meixner functions of discrete variable and 
Laguerre functions of continuous variable we substitute 7 = 0; + !, /i — 1 — h, x — s/h'm. 
the former and take the limit /i— >0,x— >oo, /ix— >s; 



M^'^)(x) 



We take also the limit of the following expressions: 
i) Differential equation 

V 1 a^-l 



(29) 



ii) Recurrence relations 



s 4 



(30) 



V(^ + « + i)(^ + i)C+i(s) 



V(^^ + a)^C-i(s) + (2n + a + 1 - s) = 0- 



(31) 



iii) Raising operator 



L+(5,n)C(s) = _i(2n + a+l-s)C(«)-«^C(«) 
= -V(n + l)(n + a; + l)C+i(s) 



(32) 
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iv) Lowering operator 



-y/n{n + a) V'n_i(s) 



(33) 



If we substitute in the differential equation (30) a = 21 + 1, u = n + I + 1, we obtain the 
reduced radial equation for the hydrogen atom, 



+ 



V 1 /(/+1) 
si s2 



u{s) = 



(34) 



the solutions of which are given by the generalized Laguerre functions 

(Z.-/-1)! 



{y^X)\ 



(35) 



This correspondence shows that we can use the difference equation of Meixner function 
as quantum model of hydrogen atom of discrete variable 



5 Calogero-Sutherland model on the lattice 

We start with the difference equation for the Hahn polynomials of discrete variable, in 
the particulr case a = /3 = A — |, namely. 



xliV-x-A-- 



A + -)(7V-i) 



(A-2,A-2^ 
li'n 



(X + 1) + 



2x (TV - A - 1) + ( A + - ) (TV - 1) 



+ n(n + 2A)/ii^ ^'^ ^\x) =0 



/il'"^''"^^(x) + 



(36) 



In the continuous limit N — > oo, — > s, the Hahn polynomials become the Jacobi poly- 
nomials of continuous variable, that in the particular case a. — = \ — \ are proportional 
to the Gegenbauer polynomials, namely. 



(37) 



The difference equation for the Hahn polynomials becomes in the continuous limit the 
differential equation for the Gegenbauer polynomials 



(38) 



Using polar coordinates s = cos q, this equation becomes. 
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with e„(A) = En{\) - Eq = [n + \f - = n(n + 2A) 

If we normalize tiie solution by the weight function p{q) = (sing)^, that is, 

V'^(g)=d„(sing)^C^(g) (39) 

we get the standard differential equation for the Calogero-Sutherland model [9] in one 
dimension 

H^^;^:{q) = K(A)V'^(g), K(A) = (n + A)" 



i/=--^-A(A-l)-V (40) 
dq^ sm^g 



6 Discrete time quantum mechanical systems 

Let H(q,p) be the time independent Hamiltonian of some quantum mechanical system in 
one dimension [3]. The Heisenberg equation of motion for the position and momentum 
operators as functions of discrete time can be written as follows 

i 1 

- (?n+l -Qn) = T; [in+l + Qn, H] (41) 

e 2 
i 1 

- iPn+l -Pn)^7, [Pn+1 + Pn, H] (42) 

e 1 

where the position operator g„ = q{ne), and momentum operator p„ = p{ne) must satisfy 
[qn,Pn] = i-yn. 

Let Pfc(x) be an polynomial of the variable x = qp + pq. It is easy to prove that 

[q,Pk{x)]^{Pk{x + 2i)-Pk{x))q 
\p,Pk{x)]^{Pk{x-2i)-Pk{x))p 

When i7 is a polynomial of the type Pk{x) we have 

i 1 

- {qn+i - qn) = [Pk{x + 2i) - Pk{x)] - (gn+i + Qn) 

e 2 

hence 

_ l-^^e {Pu{x + 2{)-P,{x)) 
l + \^e{P,{x + 2^)-P,{x)f'' 

By iteration we can calculate g„ in terms of the initial condition q^. Similarly 

" 1 + ii£ {Pk{x + 2i) - P,{x)f'' ^^^> 
Since Pk{x) is hermitian we have [qn+i,Pn+i] = [Qn,Pn]- 
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When H is & function of x we can expand it in terms of some orthornomal polynomials 
of the variable x. In particular if we take the continuous Hahn polynomials Sk{x) defined 
by the two terms recursion relation 

kSk{x) = xSk-i{x) - {k- l)Sk-2{x) 

we can express the totally symmetric polynomial Tk^kiQ^p) of all posible monomials con- 
taining k factors of q and k factors of p as [10] 

(2k)\ 



In all these cases we have 



Pn = 



l-\r:{H{:v + 2i)-H{x^ 
l + \ie {H{x + 2i) - H{x)) 

1 - ^ie {H{x + 2i) - H{x))' 



Po 



which in the limit n 



Pn 



l + ^ie {H{x + 2i) -H{x))_ 

oo, e — > 0, ne — > t, become 

— ^ q{t) = exp {itH{x)) g(0) exp {-itH{x)) 
— >■ p{t) = exp {itH (x)) p{0) exp {—itH{x)) 



(45) 
(46) 



7 Dirac equation on the lattice 

Given a function ip{nn) defined on the grid points of a Minkowski lattice with elemen- 
tary lengths e^, difference operators A(V), average operators A(V), we construct the 
Hamiltonian for the Dirac fields ipainn) on the lattice [4] 

N-l 

H = eieseg J] A1A2A3 (^m) x 

nin2n3 =0 

X 1 7o7i— A1A2A3 + 7072A1— A2A3 + 7073A1A2— A3 + moC7oAiA2A3 I V'(n^) (47) 

L ^1 ^2 ^3 J 

from which we obtain (by the Hamilton equation of motion) the Dirac equation 



i7^— A^ n ~ A0A1A2A3 ) iP{n^) = 



(48) 



the solution of which can be expressed in terms of the plane waves on the lattice, namely, 
the orthogonal functions: 



/ (^m) = n exp (-^m^nA , 



= 0, 1 . . . - 1 
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provided the dispersion relations are satisfied 



k^c = —tg 



(49) 



The transfer matrix, which carries the Dirac field from one time to the next time step can 
be obtain from the evolution operator 



Our model for the fermion field on the lattice satisfies the following conditions in order 
to escape the no-go theorem of Nielsen-Ninomiya. 

i) the Hamiltonian is traslational invariant 

ii) the Hamiltonian is Hermitian 

iii) for mo = 0, the wave equation is invariant under global chiral transformation 

iv) there is no fermion doubling 

v) the Hamiltonian is non-local (its Fourier transform has a singularity in the Brillouin 
zone) but the evolution operator, due to the Stone theorem, is unitary. 
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1 + lieoH 
1 - ^isoH 



^ (no + 1) = C/V' (no) U+ 



(50) 
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